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The IceCube Neutrino Observatory is a 1 km 3 detector currently taking data at the South Pole. 
One of the main strategies used to look for astrophysical neutrinos with IceCube is the search for a 
diffuse flux of high-energy neutrinos from unresolved sources. A hard energy spectrum of neutrinos 
from isotropically distributed astrophysical sources could manifest itself as a detectable signal that 
may be differentiated from the atmospheric neutrino background by spectral measurement. This 
analysis uses data from the IceCube detector collected in its half completed configuration which 
operated between April 2008 and May 2009 to search for a diffuse flux of astrophysical muon 
neutrinos. A total of 12,877 upward going candidate neutrino events have been selected for this 
analysis. No evidence for a diffuse flux of astrophysical muon neutrinos was found in the data set 
leading to a 90 percent C.L. upper limit on the normalization of an E~ 2 astrophysical flux of 
8.9 x 10 -9 GeV cm~ 2 s _1 sr _1 . The analysis is sensitive in the energy range between 35 TeV — 7 PeV. 
The 12,877 candidate neutrino events are consistent with atmospheric muon neutrinos measured 
from 332 GeV to 84 TeV and no evidence for a prompt component to the atmospheric neutrino 
spectrum is found. 

PACS numbers: 95.55.Vj,95.85.Ry,95.30.Cq,14.60.Lm,29.40.Ka 



I. INTRODUCTION 

There are many objects in our universe that involve 
extremely high energy processes such as matter accret- 
ing into black holes at the centers of active galaxies and 
violent explosions such as supernovae and gamma-ray 
bursts. Understanding the physics of these astrophysi- 
cal objects requires the observation of non-thermal high 
energy radiation in the form of charged cosmic rays (pro- 
tons and nuclei), gamma-rays, and neutrinos. Despite 
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progress in cosmic-ray and gamma-ray astrophysics, the 
nature of high energy astrophysical sources is still far 
from understood. Neutrinos may elucidate the funda- 
mental connection between the sources of high energy 
cosmic rays and gamma rays. 

Cosmic rays have been well studied by both space and 
ground based instruments. As astrophysical messengers, 
their main disadvantage is that they are charged particles 
and thus are deflected by magnetic fields, subsequently 
losing their directional information. High-energy gamma 
rays have been detected from many galactic and extra- 
galactic objects, but their effectiveness as cosmic mes- 
sengers over long distance scales is limited by absorption 
on extragalactic background light. Neutrinos could pro- 
vide a fundamental connection between cosmic rays and 
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gamma rays. 

Even if individual astrophysical neutrino sources are 
too weak to be detected, a superposition of all the sources 
may give rise to a detectable extraterrestrial flux. In this 
paper, we present results from a search for a diffuse flux 
of astrophysical muon neutrinos performed with the Ice- 
Cube detector using data collected in its half completed 
configuration between April 2008 and May 2009. We first 
summarize astrophysical and atmospheric neutrino mod- 
els in Section (TTJ and describe the IceCube detector in 
Section IIIII We outline in Section IIVI how our final neu- 
trino sample was obtained. The analysis methodology is 
discussed in detail in Section [V] and we present our final 
results in Section fVll 



II. ASTROPHYSICAL AND ATMOSPHERIC 
NEUTRINO FLUXES 

The benchmark diffuse astrophysical search pre- 
sented in this paper assumes an astrophysical flux, $, 
with a spectrum $ oc E~ 2 resulting from shock accelera- 
tion. In addition to the E~ 2 spectral shape, astrophysical 
models of varying normalization and spectral shapes were 
tested as well. The Waxman-Bahcall upper bound [1] was 
derived for optically thin sources assuming a $ oc E~ 2 
primary cosmic ray spectrum. Becker, Bicrmann, and 
Rhode [5] calculated the diffuse astrophysical neutrino 
flux from active galactic nuclei using observations from 
Fanaroff and Riley Class II (FR-II) radio galaxies. These 
sources were used to normalize the flux of neutrinos by 
assuming a relationship between the disk luminosity, the 
luminosity in the observed radio band, and the calculated 
neutrino flux. Mannheim [3] and Stecker [3] derived mod- 
els for optically thick Active Galactic Nuclei sources as- 
suming the sub-TeV diffuse gamma-ray flux observed by 
the Compton Gamma Ray Observatory [5] is produced 
by the decay of neutral pions. BL Lacertae objects that 
emit TeV gamma rays can be interpreted to be optically 
thin to photon-neutron interactions. The model calcu- 
lated by Miicke el al. 6\ assumes that charged cosmic 
rays are produced in these sources through the decay 
of escaping neutrons. An average spectrum of neutri- 
nos from the precursor and prompt phases of gamma-ray 
bursts is calculated in Ref. [7] by correlating the gamma- 
ray emission to the observed flux of ultra high energy 
cosmic rays. 

The primary backgrounds in the search for diffuse as- 
trophysical i/n are the atmospheric muons and neutri- 
nos arising from cosmic ray induced extensive air show- 
ers. The substantial downward-going atmospheric muon 
background persists over a wide energy range from pri- 
mary cosmic ray energies of around a GeV to the highest 
measured extensive air showers of 100 EeV [5|. These 
events were removed by using the Earth as a filter in order 
to select upward-going neutrinos traversing through the 
Earth. Two classes of atmospheric neutrinos were consid- 
ered: neutrinos arising from the decay of pions and kaons 



(the conventional atmospheric neutrino flux) and neutri- 
nos arising from the decay of charm-containing mesons 
(the prompt atmospheric neutrino flux). Detailed three- 
dimensional calculations of the energy spectrum and an- 
gular distribution of the conventional atmospheric neu- 
trino flux are summarized in Refs. [9] and |10j . The con- 
ventional atmospheric neutrino spectrum approximately 
follows an _B _3 - 7 spectrum in the TeV energy range. The 
prompt component of the atmospheric neutrino flux is 
yet to be measured, but full calculations of the prompt 
flux are given in Refs. [11H13) . The prompt component of 
the atmospheric neutrino flux is predicted to follow the 
primary cosmic ray energy spectrum which is approxi- 
mately E~ 2 - 7 . Since a hypothetical diffuse astrophysical 
neutrino flux would have a harder energy spectrum than 
atmospheric neutrino backgrounds, evidence for a diffuse 
flux would appear as a hardening of an energy-related 
observable distribution. 



III. THE ICECUBE DETECTOR 

IceCube consists of three detectors operating together. 
The main in-ice array is composed of 4800 Digital Op- 
tical Modules (DOMs) arranged in 80 strings which are 
deployed vertically with 60 DOMs per string. The detec- 
tor is deployed deep in the Antarctic ice between a depth 
of 1450 and 2450 meters. The vertical spacing between 
each DOM is 17 m and the horizontal spacing between 
each string of DOMs is 125 m giving a total instrumented 
volume of 1 km 3 . The design is optimized for the energy 
range of 100 GeV to 100 PeV [TJ]. The DeepCore exten- 
sion is deployed within the main in-ice array and consists 
of six specialized strings which lower the energy reach to 
10 GeV. IceCube was deployed in stages with the first 
string deployed during the 2005-2006 Austral summer. 
This analysis is based on one year of data taken with the 
40-string configuration (Fig. [T]) which was deployed dur- 
ing the 2007-2008 Austral summer and was operational 
from April 2008 to May 2009. 

Each DOM consists of a 13 inch (33.02 cm) pressurized 
sphere, a 10-inch (25 cm) Hamamatsu photomultiplier 
tube [T5] (model R7081-02), a mu-metal magnetic shield, 
and associated electronics responsible for the operation 
and control of the PMT as well as amplification, filter- 
ing, and calibration [TJ]. The DOMs are triggered by 
Cherenkov photons produced by charged particles in the 
Antarctic ice. In particular, ^-induced charged-current 
interactions produce muons that can traverse the entire 
IceCube array. Analog waveforms captured by the PMTs 
are digitized in situ by the DOM main board. The cap- 
ture process is initiated by a signal derived from a dis- 
criminator connected to a high-gain signal path if the 
threshold (0.25 photoelectrons) is surpassed [TB]. For the 
data set considered in this work, the triggered event was 
sent to a buffer for further filtering if it satisfied a simple 
majority trigger (SMT) of eight triggered DOMs within 
a 5 /is time window. 
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IceTop 

81 Stations, each with 

2 IceTop Cherenkov detector tanks 
2 optical sensors per tank 
324 optical sensors 



IceCube Array 

86 strings including 8 DeepCore strings 
60 optical sensors on each string 
5160 optical sensors 

December, 2010: Project completed. 86 si 



DeepCore 

strings-spaong optimized tor lower energies 
430 optical sensors 

Eiffel Tower 
324 m 



FIG. 1. Three-dimensional view of the IceCube detector lay- 
out. This work was based on the 40-string configuration which 
was half of the completed detector. Its footprint is indicated 
by the green, red, and pink circles . The 40-string configura- 
tion was operational from April 2008 to May 2009. 



Below a depth of 1450 m, the Antarctic ice is free of air 
bubbles and exhibits exceptional optical clarity with ab- 
sorption lengths ranging from 100 m to 200 m and effec- 
tive scattering lengths ranging from 20 m to 70 m [T7]H8]. 
The scattering and absorption lengths vary due to the 
concentration of dust in the glacial ice, which varies quite 
strongly with depth [19] due to varying atmospheric con- 
ditions and volcanic activity during the glacial history of 
Antarctica. The depth and wavelength dependence of the 
scattering and absorption have been measured with a va- 
riety of in-situ light sources [J7] . The ice properties have 
recently been measured over the full depth range [TS] of 
the IceCube detector using the in-situ LEDs present in 
every DOM main board resulting in what is called the 
South Pole Ice (SPICE) model. 



IV. SIMULATION AND DATA FILTERING 

A. Simulation 

This work required an accurate Monte Carlo simula- 
tion of the down-going atmospheric muon background, 
the atmospheric neutrino flux and the subsequent de- 
tector response. The simulation was used to deter- 
mine event selection criteria in order to remove the mis- 
reconstructed atmospheric muon background and in the 
profile construction method (Section [v]) to compare the 
predicted neutrino energy-correlated observable distribu- 
tion with the data to search for evidence of astrophysical 
neutrinos. 

The generation of extensive air showers initiated by 
high energy cosmic ray particles and the propagation 



of the subsequent muons through the atmosphere was 
handled by the CORSIKA (COsmic Ray Simulations for 
KAscade) JU] event generator. Hadronic interactions of 
the cosmic-ray primaries in the atmosphere were modeled 
using the SIBYLL [2V interaction model. The composi- 
tion of the primary cosmic-ray spectrum was taken from 
the Horandel poly-gonato [52] model which modeled the 
primary cosmic ray spectrum as a combination of two 
power laws for each primary particle type. 

The generation of neutrinos of all flavors was han- 
dled by the ANIS (All Neutrino Interaction Simulation) 
code [23 . ANIS uses the parton structure functions 
from CTEQ-5 [23]. Neutrinos were generated on a ran- 
dom position on the Earth's surface and then propagated 
through the Earth. The structure of the Earth is modeled 
by the PREM, or Preliminary Reference Earth Model 
[25] . In order to reduce computation time, neutrinos that 
reach the detector were forced to interact with the nearby 
Antarctic ice or bedrock to produce secondary particles 
that automatically trigger the detector. Each event was 
assigned a weight that represents the probability that this 
particular neutrino interaction occurred. Neutrinos were 
typically generated with a baseline energy spectrum of ei- 
ther E^ 1 or E~ 2 . The event weights that were calculated 
can be used to re-weight the baseline generated spectra 
to any astrophysical or atmospheric neutrino model. 

A daughter muon from a muon neutrino charged cur- 
rent interaction or an atmospheric muon passing from 
the atmosphere into earth rock was propagated using 
the Muon Monte Carlo (MMC) [25] code. MMC incor- 
porates the various continuous and stochastic energy loss 
mechanisms of ionization, bremsstrahlung, photo-nuclear 
interactions, and pair production. The Cherenkov light 
produced by the muon and the various secondaries was 
then propagated from the muon track through the detec- 
tor volume to the DOMs by using one of two methods: 
numerical tabulation and direct tracking. 

The first method was provided by the Photonics [2"T] 
software package which incorporates numerically tabu- 
lated photon distribution results of various simulation 
runs with different light sources. Photonics tables are 
computationally efficient and have the added benefit of 
allowing the full ice description to be used in the re- 
construction of muon events. The second method used 
direct photon tracking provided by the Photon Propa- 
gation Code (PPC) [TS] which allows for a more com- 
plete description of photon propagation in the Antarctic 
ice since every photon is individually tracked and propa- 
gated. This work used PPC for the simulation of neutri- 
nos and Photonics for the simulation of the background 
atmospheric muons. This choice was made since the com- 
putational efficiency of Photonics is well suited to the 
generation of a large amount of atmospheric muon back- 
ground simulation which subsequently helps to reduce 
the uncertainty of the estimated mis-reconstructed at- 
mospheric muon background. The numerical accuracy of 
PPC is appropriate for the generation of neutrino sim- 
ulation which includes the atmospheric neutrino back- 
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ground. Benchmark neutrino simulation sets generated 
with PPC and Photonics revealed that the largest dis- 
crepancy was a 30% difference in the neutrino event rate 
near a prominent dust layer 2050 m deep in the South 
Pole ice, whereas the overall neutrino event rate disagree- 
ment was 9%. 



B. Event Selection 

The event selection strategy in this analysis used the 
Earth as a filter to remove all muons from cosmic ray in- 
duced extensive air showers and retain as many neutrino- 
induced muon events as possible. The reconstructed 
energy spectrum of the neutrino sample that remained 
(roughly from the TeV to PeV energy range) was then 
analyzed using the method outlined in Section [V] for evi- 
dence of astrophysical neutrinos. The IceCube 40-string 
data set used in this analysis yielded a total live time of 
375.5 days. 

The primary trigger for this analysis was a multiplic- 
ity condition which required eight DOMs to exceed their 
discriminator threshold within a 5 /is time window. In 
addition, a local coincidence condition was enforced that 
requires the vertical neighbors of the triggered DOMs 
to trigger within 1 //s of each other. The rate for this 
primary trigger was ~ 1 kHz. Since the trigger rate 
was dominated by atmospheric muons, the data was pro- 
cessed in several stages in order to remove the atmo- 
spheric muon background and retain only neutrino events 
at the final analysis level. First, the triggered event rate 
at the South Pole was reduced to 22 Hz by using an online 
software filter. The arrival directions of the muon tracks 
in the IceCube detector were determined with a max- 
imum likelihood reconstruction procedure. The muon 
track geometry is uniquely described by an arrival direc- 
tion and a vertex position along the track which result 
in five degrees of freedom for the reconstruction. The 
likelihood function [281 parametrizes the probability of 
observing the Cherenkov photon arrival times given the 
muon track geometry. Preliminary reconstructions were 
performed using a single photoelectron (SPE) likelihood 
which utilizes the arrival time of the first Cherenkov pho- 
ton arriving in each DOM. All events reconstructed as 
upward going through the Earth (6 > 90°) were kept in 
the initial first stage of filtering. Events reconstructed 
as down-going must pass an energy cut that tightens for 
more vertical events. This ensures that truly up-going 
high energy events, initially reconstructed as down-going, 
may be correctly reconstructed and kept in the final up- 
going event sample. 

The second filtering stage involved more CPU intensive 
reconstructions performed offline outside of the South 
Pole. Among these reconstructions is the multiple photo- 
electron (MPE) fit which utilizes the likelihood descrip- 
tion of the arrival time of the first Cherenkov photon 
given N expected photons. The first photon is less scat- 
tered than the average single photon and hence the like- 



Observable and Selection Criteria 
9 > 90° 

^§1% < 8 OR ' og(L > < 7.1 

(«ch-5) (iV ch -2.5) 

a < 3° 

log(L B ayc S /L) > 25 for cos(0) < -0.2 
log(L Ba ye S /L) > (75 cos (9) + 40) for cos(0) > -0.2 
log( Z ' Bay ° 3l + LBaYCs2 ) > 35 

#SplitTime > 80° 
#SplitGeo > 80° 

NDir > 5 
LDir > 240 

|SDirj < 0.52 

TABLE I. Summary of the analysis level selection criteria 
applied to the IceCube data, neutrino simulation, and the 
atmospheric muon background simulation to obtain the final 
event sample for the analysis. 

lihood description of the detected photoelectron is mod- 
ified when this information is taken into account The 
MPE likelihood is a more sophisticated likelihood de- 
scription than the SPE likelihood reconstruction. It gives 
improved direction resolution at higher energies. Esti- 
mates of the muon energy (see the next section), the an- 
gular resolution, and quality parameters used for back- 
ground rejection are calculated during the offline process- 
ing stage. About 5% of the cosmic ray-induced muons in 
the atmosphere that trigger the IceCube detector are mis- 
reconstructed as going up through the Earth and need to 
be separated from neutrino-induced muons at the final 
analysis level. This is accomplished using quality criteria 
which are based on parameters derived from the recon- 
structed muon track. Table |T] summarizes the analysis 
cuts applied to the level 1 filtered data and simulation. 
Table |H] summarizes the number of data and simulation 
events that satisfied each successive analysis cut defined 
in Table [I] The quality parameters used to obtain the 
final analysis sample are: 

• Reconstructed zenith angle (8): The zenith an- 
gle of the reconstructed muon track is used as a cut 
parameter to select muon events with reconstructed 
directions that traverse through the Earth. 

• Reduced log-likelihood: The log-likelihood 
value of the reconstructed track was divided by 
the number of degrees of freedom of the fit. The 
number of degrees of freedom is given by the num- 
ber of triggered DOMs (N c h) minus five, which is 
the number of free parameters in the reconstruc- 
tion. Since N c h loosely correlates with the muon 
energy, the reduced log-likelihood should be ap- 
proximately energy independent. A smaller value 
indicates that the Cherenkov photons arrived at the 
individual DOMs more consistent with the likeli- 
hood description of photon arrival times. It is an 
efficient observable for separating high energy at- 
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FIG. 2. Track quality observables for data (black), atmospheric neutrino simulation (green) and mis-reconstructed atmospheric 
muon simulation (blue) after all analysis cuts have been applied. 



mospheric neutrinos from mis-reconstructed atmo- 
spheric muons. This variable was found not to be 
energy independent for lower energy atmospheric 
neutrinos, however, and was subsequently found to 
be not efficient at background rejection at lower en- 
ergies. This was resolved empirically by redefining 
the effective degrees of freedom to iV c h — 2.5 for low 
values of iV cll . 



• Error estimate from the MPE reconstruc- 
tion (a): The directional error ellipse for the MPE 
log-likelihood reconstruction was estimated follow- 
ing |29j . It provides an event by event la uncer- 
tainty of the arrival direction in the likelihood func- 
tion used in the reconstruction of muon tracks. 

• Minimum zenith angle of a two-muon re- 
construction (^spiitGco^SpiitTimc): A substantial 
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Data : 12877 Events 

Attn v H Honda 2006: 13466 Events 

Mis-reconstructed Attn. Muons : 4 Events 

Data Not Passing Bayesian Cut: 620 Events 

Atm. v, Not Passing Bayesian Cut: Events 

Atm. Muons Not Passing Bayesian Cut: 331 Events 
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FIG. 3. Reduced Log-likelihood and Bayesian likelihood ratio quality observables for data (black), atmospheric neutrino 
simulation (green), and atmospheric muon simulation (blue) after all analysis cuts have been applied. Events that do not pass 
the final two-dimensional quality criteria for the reduced log-likelihood and the Bayesian likelihood ratio are also shown. 



fraction of the atmospheric muon background re- 
sults from two or more muons triggering the Ice- 
Cube detector during the trigger window. In order 
to reduce this background, two muons were recon- 
structed for each event after splitting the triggered 
DOMs in two groups. The separation is accom- 
plished one of two ways. The first uses a geo- 
metric approach by constructing a plane perpen- 
dicular to the MPE-reconstructed track while con- 
taining the average Cherenkov photon arrival po- 
sitions. The second method is performed tempo- 
rally by using the mean Cherenkov photon arrival 
time. Each group of DOMs are used to reconstruct 
a single muon hypothesis resulting in two recon- 
structed muon tracks. Requiring the zenith angle 
from both reconstructed tracks to traverse through 
the Earth reduces the coincident atmospheric muon 
background. 

• Log-likelihood ratio between a zenith- 
weighted Bayesian reconstruction and a 
standard reconstruction: The Bayesian likeli- 
hood ratio compares the hypothesis of an up-going 



muon track with the alternative hypothesis of a 
down-going muon track consistent with the known 
zenith-dependent flux of atmospheric muons. The 
Bayesian likelihood reconstruction is performed by 
minimizing the product of the standard likelihood 
and a Bayesian prior. The Bayesian prior is based 
on the known zenith dependence of the down-going 
muon flux. Since the prior goes to zero near 
the horizon, the reconstruction always results in a 
down-going muon. Low values of the negative log- 
likelihood ratio support the alternative hypothesis 
of a down-going muon, whereas higher values indi- 
cate an up-going muon track. Further details are 
found in Ref. 28J. The likelihood ratio is zenith- 
dependent and our selection criterion based on this 
quality parameter varies with the zenith angle of 
the MPE reconstructed track. 

• Log-likelihood ratio between a zenith- 
weighted two-muon Bayesian reconstruction 
and a standard reconstruction: The two-muon 
Bayesian likelihood ratio compares the hypothesis 
of a single up-going muon track with the alter- 
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native hypothesis of two down-going muon tracks 
consistent with the known zenith-dependent flux 
of atmospheric muons. Two down-going nations 
were reconstructed separately using the DOM split- 
ting strategies discussed above. Each muon is re- 
constructed with a Bayesian prior defined with a 
zenith-dependent weight of the atmospheric muon 
flux. This observable is constructed to reject mis- 
reconstructed coincident atmospheric muons. As in 
the single muon case, low values support the alter- 
native hypothesis of two down-going atmospheric 
muons whereas higher values indicate an up-going 
muon track. 

• Number of DOMs with direct photoelec- 
trons, (NDir): The number of Cherenkov pho- 
tons arriving between —15 and +75 ns of their ex- 
pected un-scattered photon arrival times from a re- 
constructed track is known as the number of direct 
photons, or NDir [33]. More direct photons would 
indicate a better reconstructed track. 
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FIG. 4. The zenith angle distribution at final analysis level for 
375.5 days of IceCube 40-string data (black) and atmospheric 
neutrino Monte Carlo(green). 



• Direct length of the reconstructed track 
(LDir): The number of direct photons, NDir, are 
projected back onto the reconstructed track. The 
direct length, LDir |28| . is the maximum separation 
distance between these projected photons. 

• Smoothness of the reconstructed track 
(SDir): The direct photons (NDir) are again pro- 
jected back onto the reconstructed track. The 
smoothness, SDir, is a measurement of how uni- 
formly distributed these projected photons are 
along the reconstructed track. The smoothness pa- 
rameter is defined between —1 and 1. Positive val- 
ues of smoothness indicate that the projected pho- 
tons cluster at the beginning of the track, whereas 
negative values of smoothness indicate there are 
more at the end of the track. A smoothness that is 
close to indicates a uniform distribution of pro- 
jected Cherenkov photons. Further details of the 
smoothness parameter can be found in Ref . [28 . 

After all analysis level cuts have been applied, we were 
left with 12,877 candidate neutrino events below the 
horizon for the IceCube 40-string data selQ These cuts 
were designed in particular to maximize the retention ef- 
ficiency of the simulated E~ 2 astrophysical neutrino flux, 
which is 35.1% with respect to up-going events passing 
the level 1 filter. The final analysis level (after all analysis 
cuts have been applied) distributions for the track qual- 
ity observables summarized above are shown in Figs. [2][3] 
for data and Monte Carlo simulation. The zenith distri- 
bution at the final analysis level is shown in Fig. |4j The 



1 A table of the event sample used in this analysis is available at 
http: / / www.icecube.wisc.edu/science/data/ 
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FIG. 5. Effective area for + as a function of the true 
neutrino energy in intervals of the true zenith angle of the 
neutrino. The angle averaged area is represented by the solid 
black line. 



background atmospheric muon contamination was esti- 
mated to be 4 events in the final sample with a relative 
error of 60%. The background contamination was es- 
timated from simulated down-going atmospheric muons 
that survived the analysis cuts. To estimate the back- 
ground contamination, one would ideally have as much 
simulated background live-time as the data. In practice, 
the simulated background live-time was significantly less 
than the live-time of the data with eleven simulated days 
over all energies. A weighting scheme was used to in- 
crease the number of generated events at high energies 
resulting in 240 days of equivalent background live-time 
above a primary cosmic ray energy of 100 TeV. The sim- 
ulated atmospheric muons over all energy decades were 
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Purity Criterion 


Data 


Total Atm. fi 


Coincident fi 


Atm. Vy, 


E- 2 u„ 


Triggered 

bo 


3.3 x 10 10 


2.98 x 10 10 


1.72 x 10 10 


1 x 10 6 


1.03 x 10' 


LI Filter 


8.0 x 10 8 


7.5 x 10 s 


3.9 x 10 s 


1.14 x 10 5 


1,956 


9 > 90° 


2.4 x 10 8 


3.0 x 10 s 


1.79 x 10 8 


91,246 


1,353 


log(L) 


8.46 x 10 6 


4.58 x 10 6 


1.12 x 10 6 


43,183 


934 


a 


1.43 x 10 6 


1.05 x 10 6 


4.1 x 10 5 


37,174 


677 


log (Leaves /£) 
l oe ./ - L B a yc a l + - L Bay eS 2 > 


2.88 x 10 5 


2.73 x 10 5 


2.36 x 10 5 


27,411 


659 


44,309 


24,032 


17,648 


18,400 


622 


f/^i-il itTi m o 
'-'OJJ11L XIIIICJ 


22,154 


3,004 


2253 


15,771 


556 


^SplitGeo 


17,648 


1,126 


751 


15,020 


532 


NDir 


15,771 


751 


370 


14,645 


524 


LDir 


13,518 


374 


325 


14,269 


499 


SDir 


12,877 


4 





13,466 


475 



TABLE II. Number of events at each purity level for data and simulation for atmospheric muons, conventional atmospheric , 
and E~ 2 astrophysical z/ M with a normalization of N a = 10 _7 GeV cm~ 2 s _1 sr _1 for the full 40-string data set of 375.5 days. 
The background atmospheric was simulated with a total live-time of 11 days. A weighting scheme was used to increase the 
live-time at high energies resulting in 240 days of background live-time above a primary cosmic ray energy of 100 TeV. The 
quality parameter used for the purity cut is shown and the specific values of the cuts are defined in Table [I] 



then extrapolated to one year of live-time. 

The efficiency of neutrino detection for a particular 
analysis, which includes the efficiency of the analysis level 
cuts and physical effects like the absorption due to the 
Earth, can be characterized by the effective area which is 
defined as the area A e g(E, 9, 4>) of a detector that would 
have a 100% neutrino detection efficiency The total 
number of detected events is: 



iVevent, = / dE v dCL dt $ V (E V , 9, </>)A^{E, 9, <f>) (1) 



Fig. [5] shows the effective area for + as a function 
of energy for this analysis averaged over different zenith 
angle ranges. 

The final sample of candidate neutrino events is ana- 
lyzed for the presence of astrophysical neutrinos. The as- 
trophysical models considered predict a flavor ratio at the 
source of : v e : v T — 2 : 1 : 0, which subsequently os- 
cillate to a flavor ratio of : v e : v T = 1 : 1 : 1 at Earth. 
Tau neutrinos that propagate through the Earth undergo 
a regeneration effect with a branching ratio r — > fiv^Vr 
of 17% and this contribution was taken into account 
by incorporating a separate v T Monte Carlo simulation. 
The final astrophysical $ M results were derived assuming 
a flavor ratio of : v e : v T = 1 : 1 : 1 at Earth. As 
discussed in Section [TTJ, evidence for a diffuse astrophys- 
ical flux would manifest in the IceCube detector as a 
hardening at the high energy tail of the reconstructed en- 
ergy observable distribution above the expectation from 
the atmospheric spectrum. The energy-correlated ob- 
servable used in the analysis is the muon energy loss per 
unit length and is described in the next section. 



C. Energy Reconstruction 

It was natural in this analysis to use the average 
muon energy loss per meter (d_E rcco /dA") as the energy- 
correlated observable since IceCube measures the muon 
energy loss (and not the muon energy directly) in the 
form of the Cherenkov photons emitted by the vari- 
ous stochastic muon energy loss mechanisms. In order 
to estimate dE rcco /dX from the observed collection of 
Cherenkov photoelectrons (denoted by {«}) and the ex- 
pected Cherenkov photoelectron profile (denoted by fi 
and is explicitly a function of dE/dX), a log-likelihood 
based reconstruction method is used. There is an ob- 
served {n} and expected /i for every DOM in the de- 
tector. With an observed photoelectron collection {n} 
given an expected photoelectron distribution fi(dE / dX) 
binned into N bins of photoelectrons in a single DOM, a 
poisson likelihood function yields: 

log L ( d ^ c ° \{n} j = TH log Hi - Mi (2) 

where n.; and fii are the observed and expected num- 
ber of photoelectrons in the ith bin, respectively. The 
Cherenkov photoelectrons are binned according to their 
respective arrival times at the DOM. To obtain the to- 
tal likelihood function for the detector, the log-likelihood 
values of the individual DOMs were summed together: 

NdOMs 

logi t otal= ^2 logij. (3) 

i=i 

The direction and geometry of the muon were fixed 
to the results of the MPE reconstruction. The muon 
light profile fi(dE/dX) was parametrized in terms of 
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FIG. 6. The left plot shows a profile of the average reconstructed muon energy loss d_E reC o/dX vs. the true energy of the muon 
at closest approach. The error bars indicate the RMS of the reconstructed dE/dX for a given slice in the parent muon energy. 
The right plot shows a profile of the energy of the primary neutrino for different bands in reconstructed dE/dX. The error 
bars indicate the RMS of the parent neutrino energy. Shown are spectra for atmospheric neutrinos and a hypothetical E~ 2 
astrophysical z/ M flux. 




FIG. 7. Simulated neutrino energy distribution (left plot) and the simulated reconstructed muon energy loss distribution 
(right plot) of the final event sample for the Honda et. al conventional atmospheric (green) flux model, the Enberg et 
al. prompt atmospheric (light blue) flux model, and an astrophysical E~ 2 (purple) flux with a normalization of N a = 
10 _7 GeV cm" 2 s" 1 sr" 1 . 



the stochastic cascade energy which is varied until the 
likelihood function was maximized. The estimation of 
dE reco /dX is contingent on modeling /i, which depends 
both on the light yield of the muon and the optical prop- 
erties of the South Pole ice. Incorporating the muon 
light yield to the likelihood fit proved a challenge due 
to the stochastic energy loss processes of pair produc- 
tion, photo-nuclear interactions, and bremsstrahlung ra- 
diation which dominate over continuous energy losses 
above > 1 TeV. The relationship is approximately linear 
in the stochastic energy regime, however, with dE/dX = 
a + bE with a = 0.25958 GeV/mwe and b = 3.5709 x 
10~ 4 mwc -1 in ice [30] , The continuous and stochastic 
energy losses were parametrized by the coefficients a and 
b respectively and both are written in terms of meters of 



water equivalent (mwe). The reconstruction algorithm 
therefore modeled the stochastic energy loss of a muon 
as uniform along the track. 

The uniform energy loss model allows one to differen- 
tiate Eq. [2] with respect to a muon energy scale factor. 
This leads to an analytic solution for dE reco /dX in terms 
of the ratio of the total observed charge across all DOMs 
to the total predicted charge. The dE/dX reconstruc- 
tion algorithm incorporated the optical properties of the 
South Pole ice into the reconstruction. The dE/dX re- 
construction algorithm did not account for Chcrcnkov 
light from the hadronic shower initiated by the charged- 
current neutrino interaction and only reconstructed the 
energy loss due to the muon track itself. This was not 
a limitation here, since the majority of events in our fi- 
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nal sample (~ 73%) are through-going tracks where the 
neutrino interaction occurred outside the instrumented 
volume of the detector. 

The performance of the dE/dX reconstruction was 
characterized using simulated high energy muons and 
neutrinos that satisfy the analysis level selection crite- 
ria applied to the data. Since IceCube measures the en- 
ergy loss of the muon in the form of Cherenkov light 
from stochastic showers, we first characterized the in- 
trinsic resolution of the dE/dX reconstruction by using 
a sample of simulated high energy muons. The correla- 
tion of dE ICCO /dX with the muon energy closest to the 
center of the IceCube array is given in the left hand plot 
of Fig. [6] The relationship is linear over a large en- 
ergy range. The correlation changes for energies below 
1 TeV since the energy loss is no longer stochastic and 
the Cherenkov light output is nearly independent of en- 
ergy. The spread in dE TCCO /dX does not vary strongly 
as a function of the muon energy. The muon energy res- 
olution over all energy decades in the stochastic energy 
regime above a TeV of the dE/dX reconstruction is 0.27 
in log(-E). The energy resolution was estimated from the 
la width of a Gaussian fit to the reconstructed dE/dX 
distribution over all energies. The right hand plot of 
Fig. [6] shows a profile of the simulated neutrino energy 
for conventional atmospheric neutrinos and astrophysi- 
cal neutrinos vs. dE Teco /dX. Muons of a given energy 
would result in a measured dE ICCO /dX distribution with 
a mean and RMS spread as indicated by the left hand 
plot of Fig. [6] independent of the primary energy spec- 
trum. An estimate of the parent neutrino energy from 
a given measured dE rcco /dX, however, depends on the 
assumed primary energy spectrum. 

The simulated dE ICCO /dX response of the IceCube 40- 
string detector at final analysis level to the Honda et al. 
conventional atmospheric neutrino flux [10] . the Enberg 
et al. prompt atmospheric neutrino flux [11] . and a hy- 
pothetical astrophysical E~ 2 flux is shown in Fig. [7| It 
clearly shows how the different parent spectra map into 
distinguishable reconstructed muon energy spectra. 



V. ANALYSIS METHOD 

To test the compatibility of the observed dE reco /dX 
distribution in the data set with the hypotheses of muons 
arising from conventional atmospheric v^, prompt atmo- 
spheric , and astrophysical while incorporating var- 
ious sources of systematic uncertainty, we incorporated 
the frequentist approach suggested by Feldman [3T] . The 
'profile likelihood construction procedure extends the orig- 
inal frequentist method of Feldman and Cousins [32] in 
order to incorporate sources of systematic uncertainties 
parametrized as nuisance parameters in the analysis. 



A. Profile Likelihood Construction 

The profile likelihood construction method results in 
fully frequentist confidence intervals for the physics pa- 
rameters of interest (denoted by 9 T ) while using values 
of the nuisance parameters (denoted by 9 S ) that fit the 
data the best. Specifically, we first constructed a Poisson 
likelihood function and binned our dE lcco /dX observable 
distribution into N bins: 

N -n- 

L({n l }\{M,0 s )}) = U e —^ W 

where rij and /x, denote the observed and expected event 
counts in the ith dE reco /dX bin respectively. We then 
iterated over the physics parameter space and calculated 
the profile likelihood [33] ratio test statistic R p at each 
point 9 r . Defining C = — 21og(L): 

Rp(e x ) = £(9 t ,§ a ) (5) 

\L({m}\{m0 t ,§ B )}) ) 
N ( \ 

= 2 V" I fj,i - /tj + rii log I 

U v ik J 

where /2, = Hi(9 r ,8 s ) and \ii = fii(6 T ,6 s ). 9 r and 9 S de- 
note the values of the physics and nuisance parameters 
that globally minimize the profile likelihood £(6 X , 6 S ) and 
therefore describe the data the best. The value of the 
nuisance parameter 9 S conditionally minimizes the pro- 
file likelihood C(9 r , 9 S ) at the physics parameter point 8 r . 
The profile likelihood test statistic was marginalized over 
nuisance parameters in the likelihood ratio such that con- 
fidence intervals were constructed solely for the physics 
parameters of interest. 

Confidence intervals were constructed at confidence 
level a by comparing the profile likelihood test statistic 
to a critical value i? PjC rit at each point 9 r . The criti- 
cal profile likelihood value determines if a hypothesis is 
accepted or rejected at a certain confidence level. Fol- 
lowing the prescription outlined in Refs. [31] and [32] . 
we defined fl P)C rit(^r) by examining the spread of the 
profile likelihood test statistic R p (9 r ) caused by statis- 
tical fluctuations. This was facilitated by generating a 
number of Monte Carlo experiments to obtain the dis- 
tribution of i? p ((9 r ) at each physics point 9 I while fixing 

the values of the nuisance parameters to 9 S . Confidence 
intervals were constructed at confidence level a by find- 
ing the critical value of the profile likelihood, -R PlC rit(#r), 
such that the fraction a of experiments at 9 r satisfied 
R P {dr) < # P ,crit(0r)- The acceptance region is the pa- 
rameter space {9 r } such that R Pt data(0 r ) < i? P , C rit(^r) at 
a chosen confidence level a. By utilizing the profile like- 
lihood distribution to determine the confidence level, we 
have used the likelihood ratio as an ordering principle 
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in order to sort the possible experimental outcomes by 
increasing statistical significance. 



B. Systematic Uncertainties 

Systematic errors were incorporated into the analysis 
as nuisance parameters in the likelihood. The profile 
likelihood construction method discussed in the last sec- 
tion was used to define confidence regions for the physics 
parameters of interest while incorporating the various 
sources of systematic uncertainty outlined in this section. 
This work considered sources of systematic errors that af- 
fect the shape and overall normalization of the observed 
dE leco /dX distribution. The sources of systematic un- 



certainty are summarized in Table III 



One of the largest sources of systematic uncertainty 
is the overall normalization of the atmospheric neutrino 
flux. The model used in this analysis for the conven- 
tional component of the atmospheric neutrino flux was 
derived by Honda et al. |10j . where the uncertainty in 
the absolute normalization was estimated to be ±25%. 
The prompt component of the atmospheric neutrino flux 
has yet to be experimentally measured and there exists 
a large range in the overall normalization as calculated 
in various model predictions. The baseline model used 
in this analysis for the prompt component of the atmo- 
spheric neutrino flux is the calculation from Enberg et al. 
[llj where the authors quoted a standard perturbative 
QCD model prediction with an asymmetric error range 
in the overall flux normalization of —44% to +25%. 

The spectral shape of the atmospheric neutrino flux is 
affected by the uncertainty in the spectral slope of the 
primary cosmic ray spectrum. The uncertainty in the 
primary cosmic ray spectrum was estimated by consider- 
ing the uncertainty in the spectral slopes of cosmic ray 
protons (which comprise 79% of the flux) and of helium 
nuclei (15% of the flux). The remaining 6% predomi- 
nately consists of elements heavier than helium. Gaisscr 
et al. [34] estimated the spectral slope uncertainty for 
protons to be ±0.01 and for helium nuclei to be ±0.07. 
Scaling the individual spectral index uncertainties by the 
fraction of the total flux for the respective component 
gave an uncertainty in the primary cosmic ray spectral 
slope of ±0.03. This range is valid for a primary cos- 
mic ray energy up to about 1 PeV. Above this energy 
appears a region between 1 and 10 PeV known as the 
knee where the primary cosmic ray spectrum steepens 
[35j . The cosmic ray composition around the knee is still 
an active area of research. The nominal prediction of the 
model calculated in Ref. [10] was calculated up to a neu- 
trino energy of only 10 TeV. We extrapolated the Honda 
et al. model beyond this energy by assuming the model 



followed its approximate E~ 3 - 7 shape above 10 TeV. We 
note that such an extrapolation does not reflect the steep- 
ening in the spectrum expected above 1 PeV/nucleon as 
a consequence of the knee. 

There are two main sources of systematic uncertain- 
ties which affect the response of the IceCube detector 
to Cherenkov light. The first is the uncertainty in the 
absolute sensitivity of the digital optical module. The 
DOM sensitivity was assumed to be dominated by the 
uncertainty in the absolute efficiency of the photomul- 
tiplier tube which was measured to be ±8% [H)!. The 
DOM sensitivity is further reduced by a shadowing ef- 
fect from the main cable and the magnetic shield in the 
DOM, which reduces the sensitivity by 7%. The sec- 
ond dominant source of systematic uncertainty affecting 
the detector response is the uncertainty in the measured 
properties of the glacial ice at the South Pole. The mea- 
sured uncertainty in the scattering and absorption coef- 
ficients of the south pole ice was measured to be ±10% 
[18] at a flasher LED wavelength of 405 nm. 

Other relatively minor sources of systematic error were 
quantified with dedicated simulation studies. The un- 
certainty in the charged current, deep-inelastic neutrino- 
nucleon cross section was calculated in [3j5] to be ±3% 
using the parton distribution function error tables from 
[2"1] and the error calculation prescription in Ref. |37j . 
The 3% uncertainty in the cross section corresponds to 
a 3% uncertainty in the overall neutrino event rate. The 
uncertainty in the muon energy loss cross sections was 
estimated from 26 to be 1% which has a negligible ef- 
fect on the overall event rate since a decrease or increase 
in the stochastic cross sections are accompanied by a cor- 
responding increase or decrease in the muon range. The 
uncertainty in the density of the bedrock under the polar 
ice was measured to be 10% [35]. This provided a neg- 
ligible difference in the atmospheric neutrino event rates 
of < 0.1%, since the increase in the neutrino interaction 
probability is offset by a corresponding decrease in the 
range of the muon. The background contamination in 
the final event sample was estimated to be less than 1%, 
and was therefore a negligible source of systematic un- 
certainty in the analysis. 



C. Final Analysis Parameters 



The systematic errors summarized in Table III 



2 The asymmetric error range in the overall flux normalization of 
the prompt atmospheric neutrino flux only affects the overall rate 
of prompt atmospheric neutrinos 



were 

incorporated into the profile likelihood as nuisance pa- 
rameters. During minimization, each nuisance parame- 
ter was allowed to vary freely within the allowed range 
around its nominal value. The nominal values of the 
nuisance parameters correspond to the predictions of the 
Honda et al. model for the conventional atmospheric flux 
and baseline values as given by simulation for the other 
nuisance parameters. Each point in the likelihood space 
gave a specific prediction for the dE leco /dX observable 
and the profile construction method was used to define 
confidence regions for the physics parameters of interest. 
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TABLE III. Systematic uncertainties in the atmospheric neutrino flux and the detector response that affect the shape and 
overall normalization of the d-E rcC o/dX observable. The effect on the rate of triggered atmospheric is shown in the third 
column. 



Systematic Uncertainty 


Magnitude 


Atm. Rate 


Conv. Atmospheric ± Flux 


±25% 


±25% 


Prompt Atmospheric ± Flux 


-44%, +25% 


-44%, +25%0 


Cosmic Ray Spectral Slope 


±0.03 


negligible 


DOM Sensitivity 


±8% 


±15% 


Scattering &Absorption Coefficients 


±10% 


-13.5%, ±14.2% 


Neutrino-Nucleon Cross Section 


±3% 


±3% 


Muon Energy Loss 


±1% 


negligible 


Bedrock Density 


±10% 


negligible 



a The asymmetric error range in the overall flux normalization of the prompt atmospheric neutrino flux only affects the overall rate of 
prompt atmospheric neutrinos 



The likelihood methodology could be used for two main 
analyses. The primary analysis is the search for a dif- 
fuse astrophysical signal while simultaneously fitting 
for a potential prompt component of the atmospheric 
flux. In the absence of any signal, the profile likelihood 
construction could be used to measure the conventional 
atmospheric neutrino flux. 

The generic astrophysical diffuse flux was 

parametrized as an E~ 2 spectrum and the normalization 
is the main physics parameter of interest: 

$ a = N & E- 2 (6) 

where iV a has units of GeV cm~ 2 s _1 sr _1 . Astrophys- 
ical models that do not predict an E~ 2 spectral shape 
were also considered in this work. The second physics 
parameter of interest denotes the absolute normalization 
of the prompt atmospheric neutrino flux: 

/ £ \ A 7 

*P = (1 + «p) p *E„berg (7) 

\ -^median, p / 

where 1 ± a p describes the deviation in the absolute nor- 
malization of the prompt atmospheric neutrino flux from 
the reference prompt atmospheric neutrino model which 
was taken to be the calculation from Enberg et al. The 
uncertainty in the primary cosmic ray slope, A7, changes 
the shape of the predicted atmospheric neutrino flux and 
was a nuisance parameter in the analysis. It was allowed 
to float in the ±0.03 range quoted in Table [Till The shape 
dependent term was modeled by introducing an energy 
dependent weight (£ , /-E m cdian,p) A ' 7 where E medianiP is the 
median neutrino energy at final purity level. The median 
energy is 7 TeV for the nominal prediction by Enberg et 
al. 

The conventional component of the atmospheric neu- 
trino flux was treated as a nuisance parameter in the 
main analysis and was parametrized in the same fashion 
as the prompt flux: 

/ \ A 7 

*c = (1 + Oc) $Ho„da (8) 

\ -^median, c / 



where 1 ± a c describes the deviation in the absolute nor- 
malization of the conventional atmospheric neutrino flux 
from the reference model by Honda et al., A 7 is again 
the uncertainty in the primary cosmic ray slope, and 
^median, c is the median conventional atmospheric neu- 
trino energy at final purity level. A shape dependent 
term was again modeled by an energy dependent weight 
(i?/i? m0 dian,c) A7 where -E mc dian,c is 1.2 TeV for the nom- 
inal prediction by Honda et al. 

The detector efficiency, denoted by e, affects the over- 
all event rate in the IceCube detector. The magnitude 
of this systematic error combines in quadrature the sys- 
tematic uncertainties in the absolute DOM sensitivity, 
the neutrino interaction cross section, the muon energy 
loss cross sections, and the bedrock density giving an al- 
lowed range of ±8.3%. The detector efficiency nuisance 
parameter was implemented by assuming that the abso- 
lute DOM sensitivity is independent of energy. Although 
the uncertainty in the absolute DOM sensitivity affects 
the event rate for lower energy neutrino events more than 
higher energy events, this energy dependence was ne- 
glected since the primary astrophysical diffuse search is 
dominated by the high energy tail of the di? roco /dX dis- 
tribution. 

The scattering and absorption coefficients b(\ = 
405nm) and a(A = 405nm) (both measured at a LED 
wavelength of 405nm 18, 3ff]) were implemented as dis- 
crete nuisance parameters in the likelihood function. 
This was facilitated by generating Monte Carlo neutrino 
simulation sets with the scattering and absorption coef- 
ficients increased by 10%, decreased by 10%, the scatter- 
ing increased and absorption decreased by 10%, and the 
scattering decreased and absorption increased by 10%. 

To summarize, the profile likelihood used in the main 
analysis incorporated two physics parameters and five 
nuisance parameters. In the absence of any signal, 
the conventional atmospheric neutrino flux measurement 
promotes the deviation in the conventional atmospheric 
flux 1 ± a c and the uncertainty in the primary cosmic 
ray spectral slope A7 to physics parameters, giving a 
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TABLE IV. Physics parameters and nuisance parameters used 
for the astrophysical diffuse search and the measurement 
of the atmospheric neutrino spectrum 



Analysis 


Physics 


Parameters 


Nuisance Parameters 


Astro 


N a , 


1 + Q P 1 - 


|-ac, A 7 , e, 6(405), a(405) 


Atm. Vfj, 


1 + 




e, 6(405), o(405) 



likelihood with two main physics parameters and three 
nuisance parameters. The physics and nuisance parame- 
ters for the two analyses are summarized in Table [IV| 



D. Sensitivity and Discovery Potential 

A blindness procedure was followed in order to prevent 
any inadvertent tuning of the purity cuts that would bias 
the analysis. To establish a context for the unblinded re- 
sults, we quantified the limit setting potential of the anal- 
ysis (the analysis sensitivity) and the ability to discover 
an astrophysical neutrino flux (the discovery potential). 
Specifically, the sensitivity is defined as the median 90% 
upper limit obtained over an ensemble of simulated ex- 
periments with no true signal. The discovery potential 
is defined to be the strength a hypothetical astrophysical 
Vp flux required to obtain a 5er discovery in 90% of simu- 
lated experiments in the ensemble. The sensitivity of this 
analysis to a diffuse flux of astrophysical with an E~ 2 
spectrum is 1.22 x 10 -8 GeV cm~ 2 s _1 sr" 1 and the E~ 2 
discovery potential is 6.1 x 10 -8 GeV cm" 



~ 2 s 1 sr 1 . 





Data: 12877 Events 



Conv. Atmospheric v (1+o c = 0.96) 



-1.5 



-0.5 



0.5 1 1.5 

log (dE /dX)GeV/m 

3 10 reco ' 



FIG. 8. The fitted muon energy loss distribution of the final 
event sample is shown. The best fit to the data (black, shown 
with lcr error bars) consists only of conventional atmospheric 
, and no evidence is found for a prompt atmospheric flux 



or an astrophysical E 



flux. 



VI. RESULTS 




1.5 2 

N, 10* (GeV cm 2 s- 1 sr 1 ) 



After we performed the profile likelihood construction 
analysis on the d-E rcco /dA" distribution, no evidence was 
found for an astrophysical neutrino flux or a prompt 
component of the atmospheric neutrino flux. The fitted 
dE leco /dX distribution is shown in Fig. [8] and the best 
fit values of the analysis parameters to the data are sum- 
marized in Table [Vj No evidence for a signal was seen, 
so upper limits were set for astrophysical neutrino flux 
models. 



A. Upper Limits on Astrophysical Neutrino Fluxes 

The allowed regions for the normalization of the as- 
trophysical flux jV a corresponding to an E~ 2 flux 
and the normalization for prompt atmospheric neutri- 
nos are shown in Fig. [9] The upper limit for the as- 
trophysical normalization iV a at 90% confidence level 
was obtained from Fig. [9] by finding the point on the 
90% C.L. boundary along the null hypothesis of no 
prompt atmospheric neutrinos. The 90% upper limit 
on a hypothetical astrophysical <!>„ = N a E~ 2 flux at 



FIG. 9. Allowed regions for astrophysical muon neutrinos 
with an E~ 2 spectrum and prompt atmospheric neutrinos 
at 90%, 2a, and 3a confidence level. The lines indicate the 
boundary of the allowed region at the stated confidence level. 
The best fit point is indicated by the black dot at the origin. 



Earth with systematic uncertainties included is N® 0% = 
8.9 x 10~ 9 GeV cm" 2 s" 1 sr -1 . The analysis is sensitive 
in the energy range between 35 TeV and 7 PeV. The 
energy range was determined from MC simulation stud- 
ies of the analysis sensitivity, which was calculated to 
be 1.22 x 10 -8 GeV cm" 2 s" 1 sr -1 . The energy range 
was determined by introducing an energy threshold and 
ceiling such that the analysis sensitivity changed by five 
percent. The 90% upper limit derived in this work on a 
hypothetical astrophysical flux is compared to other 
limits and flux models in Fig. 
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Astrophysical neutrino models that do not predict an 
E~ 2 spectrum from various source classes were tested in 
the analysis. Of the models considered, this analysis was 
sensitive to the blazar model derived by Stecker [JS] , the 
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TABLE V. Likelihood fit results and associated errors reported by the fit. Errors are quoted as la unless otherwise noted. The 
allowed range of the nuisance parameters are also given as la Gaussian constraints. 



Parameter Best Fit 



Error 



Constrained Range 



l + a p 
1 + Qc 
A7 
c 

b c (X = 405nm) 
a(X = 405nm) 





0.96 
-0.032 

±2% 
Nominal 
Nominal 



8.9 x 10~ 9 1 eV (90% U.L.) N/A 

cm^ s sr ' ' 

0.73 (90% U.L.) N/A 

±0.16 ±0.25 

±0.014 ±0.03 

±8.3% ±8.3% 

±10% ±10% 

±10% ±10% 
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FIG. 10. Upper limits on an astrophysical flux with an E~ 2 spectrum are shown along with theoretical model predictions of 
diffuse astrophysical muon neutrinos from different sources. The astrophysical E~ 2 z/ M upper limits shown are from AMANDA- 
II [40j . ANTARES |41j . and the current work utilizing the IceCube 40-string configuration (IC40 z/ M 375.5 d). The atmospheric 
f M measurements shown are from AMANDA-II [421 143] . the IceCube 40-string (IC40) unfolding measurement [44] and the 
current work (IC40 Atmospheric Vp). 



AGN neutrino model derived by Mannheim [3], and the 
radio galaxy neutrino model from Becker, Biermann, and 
Rhode 2 . These models were rejected at the 5a confi- 
dence level. The analysis also rules out the Waxmann- 
Bahcall upper bound [I] at a 3cr confidence level. The up- 
per limits on astrophysical for the different models are 



summarized in table VI The upper limits for the models 
are expressed in terms of the model rejection factor [45j . 
which in the context of this analysis is the percentage 
of the reference model rejected at the stated confidence 
level, such that $ c .l = MRF x $ rcf . The 90% upper 
limits on these flux models are shown in Fig. [TT] We 
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TABLE VI. Upper Limits for Astrophysical f M for different Astrophysical Models. The upper limits are expressed in terms of 
the model rejection factor [45], which is the percentage of the reference model rejected at the stated confidence level such that 
$c.l = MRF x 4> rcf . 



Model 


90% C.L. 


3a C.L. 


5a C.L 


90% Energy Range (TeV-PeV) 


E- 2 ( ° eV ) 

V cm^ s sr / 


0.89 x 10" 8 


2.2 x 10" 8 


4.0 x 10~ 8 


35- 


- 7 


W-B Upper Bound 


0.4 


0.97 


1.78 


35- 


- 7 


Stecker Blazar 


0.1 


0.32 


0.42 


120 - 


- 15 


BBR FSRQ 


0.12 


0.34 


0.46 


35- 


- 7 


Mannheim AGN 


0.05 


0.21 


0.4 


9- 


1 
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FIG. 11. Upper limits derived by this work utilizing the IceCube 40-string configuration (IC40) are shown for an E~ 2 astro- 
physical flux and theoretical model predictions that do not predict an E~ 2 spectrum. 



note that the radio galaxy neutrino model from Becker, 
Biermann, and Rhode rejected at a 5<r confidence level 
was derived with a primary proton cosmic ray energy 
spectrum of E~ 2 and an optical thickness r = 0.2. The 
authors summarized in Ref. [5] a range of neutrino flux 
models with different spectral shapes of the primary pro- 
ton cosmic ray spectrum and varying optical thickness 
which are below the sensitivity of these results. 



Upper Limits on Prompt Atmospheric Neutrino 
Flux Models 



This analysis showed no evidence for a prompt com- 
ponent to the atmospheric neutrino flux. Hypotheses 
other than the reference model from Enberg et al. are 



shown in the left hand side in Fig. 12 and were tested in 
this analysis. The results of the prompt model tests are 
summarized in table |VII| and on the right hand side of 
Fig. 12 In the same fashion as the astrophysical model 
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tests described above, the upper limits on prompt atmo- 
spheric neutrinos were expressed in terms of the model 
rejection factor. The standard calculation from Enberg 
et al. which is used as the reference flux in this analysis 
was rejected at 90% confidence level valid from an energy 
range between 9 TeV and 613 TeV. 

Under the assumptions of the present analysis, we 
reject the RQPM [12] (Recombination Quark Parton 
Model) at a 3er confidence level, which strongly disfa- 
vors the authors' non perturbative quantum chromody- 
namics (pQCD) approach to calculating the prompt flux. 
The maximum and minimum calculations from Enberg 
et al. represent an allowed theoretical uncertainty band 
due to the authors' choices of the parton distribution 
function (PDF), the charm quark mass, and the factor- 
ization scale which affect the pQCD calculation of the 
prompt atmospheric neutrino flux. The reference model 
used the MRST 2001 g7] for the PDF, a factorization 
scale /iF = 2m c where m c is the charm quark mass, and 
a charm quark mass of 1.3 GeV. The theoretical uncer- 
tainty represented by the minimum and maximum calcu- 
lations were obtained by varying the quark mass between 
1.3 GeV and 1.5 GeV, varying the factorization scale be- 
tween fip = m c and fip = 2m c , and varying the PDFs 
by using MRST 2001 or CTEQ 6 0S]. We ruled out 
the maximum calculation at 95% C.L. and the standard 
prediction at 90% C.L. These limits favor the CTEQ 6 
parameterization of the PDF, a lower quark mass, and a 
low factorization scale. 

TABLE VII. Upper limits on prompt atmospheric neutrinos 
for different models. The upper limits are expressed in terms 
of the model rejection factor [45], which is the percentage 
of the reference model rejected at the stated confidence level 
such that $c.l = MRF x $ re f . 



Model 


90% 


95% 


3(7 


90% Energy 


Range 


Enberg (Minimum) 


1.25 


1.8 


3.6 


9- 


615 


Enberg (Standard) 


0.73 


1.1 


2.2 


9- 


613 


Enberg (Maximum) 


0.53 


0.85 


1.89 


9- 


610 


Naumov RQPM 


0.2 


0.41 


0.87 


9- 


620 


Martin MRS 


2.1 


4.0 


8.9 


9- 


613 



(TeV) 



C. Measurement of the Atmospheric Neutrino 
Spectrum 

There was no evidence for astrophysical neutrinos or 
prompt atmospheric neutrinos in the final event sample, 
and therefore the final neutrino distribution was inter- 
preted as a flux of conventional atmospheric muon neu- 
trinos. The profile construction method was used to mea- 
sure the atmospheric neutrino flux in order to determine 
the normalization and any change in shape from the ref- 
erence atmospheric neutrino flux model considered. The 
best fit result of the atmospheric neutrino flux is of the 



form: 

/ P \ -0.032±0.014 
<&Bc S tFit = (0.96 ± 0.16) ( 1 - J $Honda (9) 

where the normalization of the atmospheric neutrino flux was 
found to be 4% ±16% lower than the nominal prediction from 
Honda et al. and the spectral index was found to be steeper 
by A7 = -0.032 ± 0.014. The allowed regions of (1 + a c ) and 
A7 are shown with the band of allowed atmospheric neutrino 
spectra in Fig. |13| The displayed band of allowed atmospheric 
neutrino spectra in black was constructed from the envelope of 
the set of curves allowed by the 90% boundary in the left plot 
of Fig. |13| We note that our best fit for the spectral index 
rejected the nominal spectral index prediction from Honda 
et al. at a 95% confidence level, neglecting the theoretical 
uncertainty. The overall flux normalization is consistent with 
the nominal model prediction at 90% confidence level. The 
energy range of the atmospheric neutrino flux measurement is 
valid from 332 GeV to 84 TeV. This energy range was derived 
from the median atmospheric neutrino energy as predicted by 
simulation for the lowest and highest dE leco /dX values from 
the data. 

Also shown in the right plot of Fig. [13] is the atmospheric 
neutrino unfolding analysis discussed in Ref. [44] where no 
prior assumption was made regarding the shape of the atmo- 
spheric neutrino spectrum, while this work parametrized the 
atmospheric neutrino flux as a power law. Such a bias in the 
prior assumption of the atmospheric neutrino flux resulted in 
tighter error bars than in the unfolding analysis. 



VII. CONCLUSIONS AND OUTLOOK 

To summarize, we have set the field's most stringent limit 
on astrophysical muon neutrinos from unresolved sources. 
The 90% upper limit on an astrophysical flux with an E~ 2 
spectrum is 8.9 x 10~ 9 GeV cm -2 s _1 sr _1 , valid from the 
energy range of 35 TeV to 7 PeV. Several optimistic astro- 
physical neutrino production models have been rejected at a 
5a confidence level. We have set limits on the prompt com- 
ponent of the atmospheric neutrino flux, with a preference for 
perturbative QCD models in the energy range between 9 TeV 
and 613 TeV. Finally, we have also measured the atmospheric 
muon neutrino flux from 332 GeV to 84 TeV and found a 
fit to the overall normalization of the atmospheric neutrino 
flux that is 4% lower than the calculation from Honda et al. 
[10] . The preferred spectrum is somewhat steeper than the 
assumed extrapolation of the Honda spectrum, perhaps re- 
flecting the steepening of the spectrum associated with the 
knee, as discussed in Ref. Overall, our result here is 

consistent with other measurements made of the atmospheric 
neutrino flux with the IceCube detector. The 90% error band 
on the measured flux overlaps with the 90% error band of the 
result from IceCube's predecessor, the AMANDA-II experi- 
ment |4"2"] . 

The 90% upper limits on the astrophysical i/ M models and 
the prompt component of the atmospheric neutrino flux are 
dependent on the assumptions made on the conventional at- 
mospheric neutrino flux. We extrapolated the model from 
Honda et al. above the maximum calculated energy of 10 TeV 
by assuming the conventional spectrum continued to follow its 
approximate E~ 3,7 spectral shape. This extrapolation is de- 
pendent on the location and primary cosmic-ray composition 
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FIG. 12. The left plot shows the predicted prompt atmospheric neutrino fluxes averaged over zenith angle and multiplied by 
E :i to enhance features. The Honda 2006 model is shown for comparison. The right plot shows the 90% confidence level upper 
limit on the prompt models. 




FIG. 13. The left plot shows allowed regions for the normalization (1 + a c ) and the change in spectral index (A7) of the 
conventional atmospheric neutrino flux relative to Honda et al [10) . The right plot compares the angle-averaged + 
measurement of the atmospheric neutrino flux of this work to the model prediction from Honda et al. The fluxes are multiplied 
by E 3 to enhance features. The displayed set of black curves is the band of allowed atmospheric neutrino spectra constructed 
from the 90% boundary of the left plot in this work as described in the text. Also shown (blue triangles) is the IceCube 
40-string unfolding analysis |44j . and the AMANDA-II (empty brown band) result [42] . 



of the knee which are currently not well known. The spec- 
trum of all nuclei steepens in an energy region around 3 PeV 
total energy per nucleus. With standard assumptions about 
the composition, and assuming that the underlying physics 
depends on magnetic rigidity, the spectrum of protons must 
become steeper around 1 PeV or lower. This has the conse- 
quence that the spectrum of neutrinos from decay of pions 
and kaons must steepen at around 100 TeV [49], which is 
just in the crossover region for the prompt component. What 
is needed is a detailed calculation of atmospheric neutrinos 
based on a realistic treatment of the primary cosmic-ray spec- 
trum and composition through the knee region that extends 
beyond the current limits of 10 TeV. Because of the steepen- 
ing, the limits on prompt neutrinos in Fig. [12] will be relaxed 
to some extent. For future work, it is critically important also 
to obtain more precise measurements of the primary compo- 



sition and spectra in the knee region. The KASCADE ex- 
periment [50] already suggests that the proton component is 
suppressed in the knee region. IceCube itself, with its surface 
component IceTop, has the potential to measure the spectrum 
and composition in the knee region and beyond. 

The stringent 90% upper limit on a diffuse astrophysical 
flux of muon neutrinos reported by this work implies that the 
IceCube detector in its 40 string configuration (as used in this 
analysis) is not yet sensitive enough to discover astrophysical 
neutrinos from unresolved sources. The full 86-string array 
has been completed during the 2010-2011 summer construc- 
tion season at the South Pole. A 5a discovery of an astro- 
physical E~ 2 Vfj, flux at the 90% limit derived by this work 
will take three years of the full IceCube array. 

This time scale for discovery can be made shorter by an 
improved understanding of the various sources of systematic 
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uncertainty and considering new analysis techniques. With 
a proper measurement of the prompt component of the at- 
mospheric neutrino flux, the time scale for discovery becomes 
more tractable. In this analysis we have not yet made use of 
the difference in angular behavior of prompt neutrinos (which 
are isotropic) and conventional atmospheric neutrinos (which 
have a higher intensity near the horizon) . Analyses dedicated 
to the study of leptons from the decay of charmed mesons 
would also yield a better understanding of the physics of air 
showers and atmospheric neutrinos. Strategies other than 
using atmospheric to search for the prompt component in- 
volve a thorough investigation of the down-going muon flux 
and a measurement of the atmospheric neutrino spectrum 
from v c . The measurement of the atmospheric v c flux has 
an advantage that the transition energy from conventional v e 
to prompt v c occurs at an order of magnitude lower in energy 
than in v^. 

The event selection in this analysis used the Earth as a fil- 
ter to remove the large down-going atmospheric muon back- 
ground. An improved simulation of atmospheric muons would 
allow a diffuse astrophysical search to incorporate the 
down-going region in the analysis and search for astrophysical 
neutrinos over the entire sky. We note that there is a slight tilt 
in the measured angular distribution of atmospheric neutrinos 
with respect to our extrapolation based on the angular depen- 
dence in [10], as illustrated in Fig. [4] This discrepancy did 
not affect our limit on astrophysical or our reconstructed 
atmospheric neutrino spectrum. However, understanding the 
origin of the discrepancy is important for future work. 

Although this analysis focuses on v^, IceCube is sensitive 
to all flavors of neutrinos. As the detector grows, reconstruc- 
tion methods mature, and the understanding of the various 
sources of systematic uncertainty improve, it would be possi- 
ble to combine event topologies from different neutrino flavors 



in a multi-flavor analysis. A simultaneous search for neutri- 
nos of all flavors from unresolved astrophysical sources would 
be significantly more sensitive than an analysis focusing ex- 
clusively on a single neutrino flavor. 
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